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Subwavelength aperture arrays in thin metal films can enable enhanced transmission of light and 
matter (atom) waves. The phenomenon relies on resonant excitation and interference of the plasmon 
or matter waves on the metal surface. We show a new mechanism that could provide a great resonant 
and nonresonant transmission enhancement of the light or de Broglie particle waves passed through 
the apertures not by the surface waves, but by the constructive interference of diffracted waves 
(beams generated by the apertures) at the detector placed in the far-field zone. In contrast to other 
models, the mechanism depends neither on the nature (light or matter) of the beams (continuous 
waves or pulses) nor on material and shape of the multiple-beam source (arrays of 1-D and 2-D 
subwavelength apertures, fibers, dipoles or atoms). The Wood anomalies in transmission spectra of 
gratings, a long standing problem in optics, follow naturally from the interference properties of our 
model. The new point is the prediction of the Wood anomaly in a classical Young-type two-source 
system. The new mechanism could be interpreted as a non-quantum analog of the superradiance 
emission of a subwavelength ensemble of atoms (the light power and energy scales as the number of 
light-sources squared, regardless of periodicity) predicted by the well-known Dicke quantum model. 
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The scattering of waves by apertures is one of the basic 
phenomena in the wave physics. The most remarkable 
feature of the light scattering by subwavelength aper- 
tures in a metal screen is enhancement of the light by 
excitation of plasmons in the metal. Since the obser- 
vation of enhanced transmission of light through a 2D 
array of subwavelength metal nanoholes the phe- 
nomenon attracts increasing interest of researchers be- 
cause of its potential for applications in nanooptics and 
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enhancement of light is a process that can include res- 
onant excitation and interference of surface plasmons 
SHU, Fabry-Perot-like intraslit modes @, 0, E S Gil , 
and evanescent electromagnetic waves at the metal sur- 
face In the case of a thin screen whose thickness 
is too small to support the intraslit resonance, the ex- 
traordinary transmission is caused by the excitation and 
interference of plasmons on the metal surface [1, 0, @. 
Recently, the enhanced transmission through subwave- 
length apertures by excitation of the matter- wave analog 
of surface plasmons was predicted also for de Broglie mat- 
ter waves |l2f . For some experimental conditions, many 
studies [IMJ, M, 0, M, H 0, EH, [H 0, [H] in- 
dicated an essential role of the surface plasmons in the 
enhancement of light waves. For an example, the study 
[1^ showed that a perfect conductor whose surface is 
patterned by an array of holes can support surface po- 
laritons, which just mimic a surface plasmon. Nowadays, 
it is generally accepted [25| that the excitation and in- 
terference of surface plasmon-polaritons play a key role 
in the process of enhancement of light waves in the most 
of experiments (also, see a recent comprehensive reviews 
2(| 27 1 ) . In the present study [28| , we show a new mecha- 



nism that could provide a great resonant and nonresonant 
transmission enhancement of the light or de Broglie parti- 
cle waves passed through the apertures not by the surface 
waves, but by the constructive interference of diffracted 
waves (beams generated by the apertures) at the detector 
placed in the far-field zone. 

The transmission enhancement by the constructive in- 
terference of diffracted waves at the detector can be ex- 
plained in terms of the following theoretical formulation. 
We first consider the transmission of light through a 
structure that is similar, but simpler than an array of 
holes, namely an array of parallel subwavelength-width 
slits in the metal screen. In some respects, the resonance 
excitation and interference of surface plasmon-polaritons 
in these two systems are different from each other [29j. 
The difference, however, is irrelevant from a point of 
view of our model. Indeed, the excitation of plasmon- 
polaritons and coupling between the apertures do not 
affect the principle of the enhancement based on the con- 
structive interference of diffracted waves (beams gener- 
ated by the independent apertures) at the detector placed 
in the far-field zone. The resonant excitation of the plas- 
mons or trapped electromagnetic modes, as well as the 
coupling between apertures could provide just additional, 
in comparison to our model, enhancement by increasing 
the power (energy) of each beam. Therefore, our model 
considers an array of slits, which are completely inde- 
pendent from each other. We also assume, for the sake 
of simplicity, that the metal is a perfect conductor. Such 
a metal is described by the classic Drude model for which 
the plasmon frequency tends towards infinity. The beam 
produced by each independent slit is found by using the 
Neerhoff and Mur model, which uses a Green's function 
formalism for a rigorous numerical solution of Maxwell's 
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equations for a single, isolated slit 

Em mm sim. 

In the model, the screen placed in vacuum is illuminated 
by a normally incident TM-polarized wave with the wave- 
length A = 2itc/uj = 2-jr/k. The magnetic field of the 
incident wave H(x,y, z,t) — U(x)exp(— i(kz + u>t))e v is 
supposed to be time harmonic and constant in the y 
direction. The transmission of the slit array is deter- 
mined by calculating all the light power of the ensem- 
ble of beams in the observation plane. To clarify the 
numerical results, we then present an analytical model, 
which quantitatively explains the resonant and nonreso- 
nant enhancement in the intuitively transparent terms of 
the constructive interference of diffracted waves (beams 
generated by the apertures) at the detector placed in 
the far-field zone. Finally, we show that the mechanism 
depends neither on the nature (light or matter) of the 
beams (continuous waves or pulses) nor on material and 
shape of the multiple-beam source (arrays of 1-D and 2-D 
subwavelength apertures, fibers, dipoles or atoms). 

Let us first investigate the light transmission versus 
the wavelength by using the rigorous numerical model. 
The model considers an ensemble of M waves (beams) 
produced by M independent slits of width 2a and period 
A in a screen of thickness b. The transmission of the 
slit array is determined by calculating all the light power 
P(A) radiated by the slits into the far-field diffraction 
zone, oo, co] at the distance z>A from the screen. 
The total per-slit transmission coefficient, which repre- 
sents the per-slit enhancement in transmission achieved 
by taking a single, isolated slit (beam) and placing it in 
an M-slit (M-beam) array, is then found by using an 
equation Tm(A) = P{X)/MP\, where P\ is the power 
radiated by a single slit. Figure 1 shows the transmis- 
sion coefficient Tm(A), in the spectral region 500-2000 
nm, calculated for the array parameters: a = 100 nm, 
A = 1800 nm, and b = 5 x 10~ 3 \ max . The transmit- 
ted power was computed by integrating the total energy 
flux at the distance z = 1 mm over the detector region 
of width Aa; = 20 mm. The transmission spectra Tm(A) 
is shown for different values of M. We notice that the 
spectra Tm(A) is periodically modulated, as a function of 
wavelength, below and above a level defined by the trans- 
mission Ti(A) = 1 of one isolated slit. As M is increased 
from 2 to 10, the visibility of the modulation fringes in- 
creases approximately from 0.2 to 0.7. The transmission 
Tm exhibits the Fabry-Perot like maxima around wave- 
lengths A„ = A/n. The spectral peaks increase with 
increasing the number of slits and reach a saturation 
{jrnax _ 5 ) in amp ii tu de by M = 300, at A w 1800 
nm. The peak widths and the spectral shifts of the reso- 
nances from the Fabry-Perot wavelengths decrease with 
increasing the number M of beams (slits). An analy- 
sis of Fig. 1, indicates that the power (energy) enhance- 
ment and dispersion are the general interference proper- 
ties of the ensemble of beams. Therefore, the enhance- 
ment and suppression in the transmission spectra could 
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FIG. 1: The per-slit transmission TW(A) of an array of in- 
dependent slits of period A in a thin (b <C A) screen versus 
the wavelength for the different number M of slits. There 
are three Fabry-Perot like resonances at the wavelengths 
AnPsA/n, where n=l, 2 and 3. 



be considered as the natural properties also of the peri- 
odic array of independent subwavelength slits. The spec- 
tral peaks are characterized by asymmetric Fano-like pro- 
files. Such modulations in the transmission spectra are 
known as Wood's anomalies. The minima and maxima 
correspond to Rayleigh anomalies and Fano resonances, 
respectively (3rj| . The Wood anomalies in transmission 
spectra of gratings, a long standing problem in optics, 
follow naturally from the interference properties of our 
model. The new point, in comparison to other mod- 
els 



371 l38l|). is the prediction of a weak Wood anomaly 



in a classical Young-type two-source system (see, Fig. 1). 

The above-presented analysis is based on calculation of 
the energy flux of a beam array, in which the electromag- 
netic field of a single beam is evaluated numerically. The 
transmission enhancement and dispersion were achieved 
by taking a single, isolated slit (beam) and placing it in 
a slit (beam) array. The interference of diffracted waves 
(beams generated by the slits) at the detector placed 
in the far-field zone could be considered as a physical 
mechanism responsible for the enhancement and disper- 
sion. To clarify the results of the computer code and 
gain physical insight into the enhancement mechanism, 
we have developed an analytical model, which yields sim- 
ple formulas for the electromagnetic field of the beam 
produced by a single slit. For the field diffracted by a 
narrow (2a <C A, b > 0) slit into the region \z\ > 2a, 
the Neerhoff and Mur model simplifies to an analytical 
one [3§|. For the magnetic H — (0,H y ,0) and electric 
E = (E x , 0, E z ) components of the single beam we found 
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the following analytical expressions: 

H v (x, z) = iaDF^(k[x 2 + z 2 } 1/2 ), 



(1) 



E x (x, z) = ~az[x 2 + z 2 ]-^ 2 DFl(k[x 2 + z 2 } 1 / 2 ), (2) 



and 

E z (x,z) = ax[x 2 + z^^DFlmx 2 + z 2 ] 1/2 ), 
where 

D = 4fc- 1 [[exp(ifc6)(aA - fc)] 2 - (aA + k) 2 }- 1 

and 



(3) 



(4) 



A = F^(ka) + -[F Q {ka)Fl{ka) + F 1 (ka)F^(ka)]. (5) 

Here, Fq, F and F\ are the Hankel and Struvc 
functions, respectively. The beam is spatially inhomo- 
geneous, in contrast to a common opinion that a sub- 
wavelength aperture diffracts light in all directions uni- 
formly 40]. The electrical and magnetic components of 
the field produced by a periodic array of M independent 
slits (beams) is given by E(x, z) = ^ m=1 E m (x + mA, z) 
and H (x, z) = 2m=i H m (x + mA, z), where E rn and H m 
are the electrical and magnetic components of the ru- 
th beam generated by the respective slit. As an exam- 
ple, Fig. 2(a) compares the far-field distributions E and 
H calculated by the analytical formulae (1-5) to that 
obtained by the rigorous computer model. We notice 
that the distributions are undistinguishable. The field 
power P{E, H) is found by integrating the energy flux 
S = E x H* + E* x H. Therefore, the analytical model 
accurately describes also the coefficient Tm of the system 
of M independent subwavelength slits (beams). The an- 
alytical model not only supports results of our rigorous 
computer code (Fig. 1), but presents an intuitively trans- 
parent explanation (physical mechanism) of the enhance- 
ment and suppression in transmission spectra in terms of 
the constructive or destructive interference of the waves 
(beams produced by the subwavelength-width sources) 
at the detector placed in the far-field zone. The array- 
induced decrease of the central beam divergence by the 
far-field multiple-beam interference (Figs. 2(b) and 2(c)) 
is relevant to the beaming light [4ll | , as well as the non- 
diffractive light and matter beams [42j, [43j . The ampli- 
tude of a beam (evanescent spherical-like wave) produced 
by a single slit rapidly decreases with increasing the dis- 
tance from the slit (1-3). However, due to the multiple- 
beam interference mechanism of the enhancement and 
beaming, the array produces in the far-field zone a prop- 
agating wave with low divergence. Such a behavior is 
in agreement with the Huygens-Fresnel principle, which 
considers a propagating wave as a superposition of sec- 
ondary spherical waves. 
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FIG. 2: The electrical and magnetic components of the field 
produced by an array of M independent slits (beams), (a) 
The distributions Re(E x (x)) (A and D), Re(H y (x)) (B and 
E), and Re(10£ z (as)) (C and F) calculated for M = 10 and A 
= 1600 nm. The curves A, B, and C: rigorous computer code; 
curves D, E, and F: analytical model, (b) Ke(E x (x)) for 
M=l: analytical model, (c) Re^^a;)) for M=5: analytical 
model. 



We now consider the predictions of our analytical 
model in light of the key observations published in the lit- 
erature for the two fundamental systems of wave optics, 
the one-slit and two-slit systems. The major features of 
the transmission through a single subwavelength slit are 
the intraslit resonances and the spectral shifts of the res- 
onances from the Fabry-Perot wavelengths Q ■ In agree- 
ment with the predictions [3] , the formula (4) shows that 
the transmission T = P/P = (a/k)[Re(D)} 2 + [Im(D)] 2 
exhibits Fabry-Perot like maxima around wavelengths 
A rl = 2b/ n, where Pq is the power impinging on the 
slit opening. The enhancement and spectral shifts are 
explained by the wavelength dependent terms in the de- 
nominator of Eq. (4). The enhancement (T(Ai)«6/7ra 
[39l ]) is in contrast to the attenuation predicted by the 
model Q. Although, our model considers a screen of 
perfect conductivity, polarization charges develop on the 
metal surface. The surface polaritons do not adhere 
strictly to traditional surface plasmons. Nevertheless, 
at the resonant conditions, the system redistributes the 
electromagnetic energy by the surface polaritons in the 
intra-slit region and around the screen. Thus, additional 
energy could be channeled thought the slit in compar- 
ison to the energy impinging on the slit opening. The 
mechanism is somewhat similar to that described in the 
study 19]. This study showed that a perfect conduc- 
tor whose surface is patterned by an array of holes can 
support surface polaritons that mimic a surface plas- 
mon in the process of channeling of additional energy 
into the slit. We considered TM-polarized modes be- 
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FIG. 3: The per-slit transmission coefficient T(X) versus 
wavelength for the Young type two-slit experiment [371 ]. Solid 
curve: experiment; dashed curve: analytical model. Parame- 
ters: a = 100 nm, A = 4900 nm, and b = 210 nm. 



cause TE modes are cut off by a thick slit. In the case 
of a thin screen, TE modes propagate into slit so that 
magneto-polaritons develop. Because of the symmetry 
of Maxwell's equations the scattering intensity is for- 
mally identical with E and H swapping roles. Again, 
the magneto-polaritons could provide channeling of ad- 
ditional energy into the slit. This enhancement mech- 
anism is different from those baced on the construc- 
tive interference of the waves (beams produced by the 
subwavelength-width sources) at the detector placed in 
the far-field zone. The Young type two-slit (two-beam) 
configuration is characterized by a sinusoidal modulation 
of the transmission spectra (for an example, see 12(A) 



Refs. 37, \3 



The modulation period is inversely 
proportional to the slit separation A. The visibility V 
of the fringes is of order 0.2, independently of the slit 
separation. In our model, the transmission T2 depends 
on the interference-like cross term J'[F 1 1 (a:i)[iF 1 (x2)]* + 
Fl{xi)*iF^(x2)]dx, where x\ = x and £2 = x + A. 
The high-frequency interference-like modulations with 
the sideband-frequency /«(A) ~/i(A) + / 2 (A, A)~l/A 
(Figs. 1 and 3) are produced like that in a classic hetero- 
dyne system by mixing two waves having different spatial 
frequencies, /1 and ji- Although our model ignores the 
enhancement by the plasmon-polaritons, its prediction 
for the transmission (T 2 mox «l.l), the visibility (V&0.1) 
of the fringes and the resonant wavelengths A n « A/n 
compare well with the plasmon-assisted Young's type ex- 
periment [37} (Fig. 3). It should be noted that in the 
case of b > A/2, the far- field interference resonances at 
A„ w A/n could be accompanied by the intraslit polari- 
ton resonances at A„ w 2b /n. One can easily demon- 
strate such behavior by using the analytical formulas (1- 
5). The interference of two beams at the detector is not 
only the contribution to enhanced transmission. There 
could be enhancement also due to the energy redistribu- 
tion by the resonant intraslit plasmon-polaritons and/or 
by the surface waves with resonant coupling through the 
slits. We stress, however, that the plasmons or trapped 
electromagnetic modes do not affect the principle of the 
enhancement based on the constructive interference of 
diffracted waves (beams generated by the independent 
subwavelength-width apertures) at the detector placed 
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FIG. 4: The per-slit transmission Tm(A) versus wavelength 
for the different values of A and M: (A) A = 100 nm, M = 2; 
(B) A = 500 nm, M = 2; (C) A = 3000 nm, M = 2; (D) A = 
100 nm, M = 5; (E) A = 500 nm, M = 5; (F) A = 3000 nm, 
M = 5. Parameters: a — 100 nm and b = 10 nm. There are 
two enhancement regimes at A < A and A>A. 



in the far-field zone. The plasmon-polaritons could pro- 
vide just additional enhancement by increasing the power 
(energy) of each beam. This kind of enhancement is 
of different nature compared to our model, because the 
model requires nether resonant excitation of the intraslit 
plasmon-polaritons nor coupling between the slits (see, 
also Refs. 0, 0, S3). 

In order to gain physical insight into the mechanism of 
plasmonless and polaritonless enhancement in a multiple- 
slit or multiple-beam (M>2) system, we now consider the 
dependence of the transmission Tm(A) on the slit (beam) 
separation A. According to the Van Citter-Zernike coher- 
ence theorem, a light source (even incoherent) of radius 
r = M (a + A) produces a transversally coherent wave 
at the distance z<TtRr/\ in the region of radius R. In 
the case of A <C A, the collective coherent emission of an 
ensemble of slits (beams) generates the coherent electro- 
magnetic field {E — ^2 m=1 E m exp(i(p m )^sM Eiexp(itp) 
and H&MHiexp(i(p)) in the far-field zone of the region 
of radius R = 00. This means that the beams arrive to 
the detector with the nearly same phases <f m {x) ~ 
(see, also Ref. [HI). Consequently, the beams add coher- 
ently and the power (energy) of the emitted light scales 
as the number of beams squared, regardless of period- 
icity, P w M 2 P\. Thus, the transmission enhancement 
(Tm = P/MPx) grows linearly with the number of slits, 
Tm~M. For a given value of M, in the case of A < A, 
the transmission Tm(A) monotonically (non-resonantly) 
varies with A (see, Fig. 4). At the appropriate condi- 
tions, the transmission can reach the 1000-times non- 
resonant enhancement (M = \z/-KR(a + A)). In the 
case of R > Xz/nr or A>A, the beams arrive to the 
detector with different phases (p m (x). Consequently, the 
power and transmission enhancement grow slowly with 
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the number of beams (Figs. 1-4). The constructive or 
destructive interference of the beams leads respectively 
to the enhancement or suppression of the transmission 
amplitudes. Although, the addition of beams is not 
so efficient, the multiple-beam interference leads to en- 
hancements and resonances (versus wavelength) in the 
total power transmitted. In such a case, the transmis- 
sion coefficient Tm exhibits the Fabry-Perot like maxima 
around the wavelengths A„ = A/n. We stress again that 
the constructive or destructive interference of beams at 
the detector requires nether the resonant excitation of 
plasmon-polaritons nor the coupling between radiation 
phases of the slits. The plasmon-polariton effects could 
provide just additional enhancement by increasing the 
power (energy) of each beam. Our consideration of the 
subwavelength gratings is similar in spirit to the dynam- 
ical diffraction models [13|, Airy- like model I14II . and es- 
pecially to a surface evanescent wave model [llj. In the 
case of A >> A, our model is in agreement with the the- 
ories of conventional (non-subwavelength) gratings |47j |. 

In the above-presented multiple-beam interference 
model, we have considered a particular light-source, 
namely an array of subwavelength metal slits. One can 
easily demonstrate the interference mediated enhance- 
ment and suppression in the transmission and reflec- 
tion spectra of an arbitrary array of subwavelength- 
dimension sources of light or de Broglie particle waves 
by taking into account the interference properties of 
Young's double-source system. At the risk of belabor- 
ing the obvious, we now describe the phenomenon. In 
the far-field diffraction zone, the radiation from two pin- 
holes of Young's setup is described by two spherical 
waves. The light intensity at the detector is given by 
1(f) = \(E/ri)exp(ikri + ipi) + (E/r 2 )exp(ikr 2 + ip 2 )\ 2 = 
h + 1 2 + 2(I 1 I 2 ) 1 / 2 cos([kr 1 + (fn] - [kr 2 + tp 2 }). The corre- 
sponding energy isW = J J[Ii+I 2 + 2(IiI 2 ) 1 ^ 2 cos([kri + 
(pi] — [kr 2 + ip 2 ])]dxdy. Here, we use the units At = 1. 
In conventional Young's setup, which contains the pin- 
holes separated by the distance A >> A, the interference 
cross term (energy) vanishes. Therefore, the energy is 
given by W = J J{h + h)dxdy = W x +W 2 = 2W , 
where W\ — W 2 — Wo- In the case of Young's sub- 
wavelength system (A << A, correspondingly r\ = r 2 
for any coordinate x or y), the energy W = W\ + W 2 + 
2 J r(/ 1 / 2 ) 1 / 2 cos((/5i — <p 2 )dxdy. The first-order corre- 
lation term could provide the enhancement or suppres- 
sion of both the intensity and energy of the light field 
at the detector (see, also Ref. 45j, |48(). Indeed, at the 
phase condition tfi — tp 2 — 0, the energy enhancement 
is given by W — 4Wo- In the case of ip\ — ip 2 = n, 
the destructive interference of the two waves leads to the 
zero transmission, W = 0. The same phase conditions 
provide the enhancement or suppression of transmitted 
energy by quantum two-source interference (for exam- 
ple, see formulas 4. A. 1-4. A. 9 [49]). The enhancement or 
supression by the classic or quantum interferece at the 



detector depends neither on the nature (light or matter) 
of the beams (continuous waves or pulses) nor on ma- 
terial and shape of the multiple-beam source (arrays of 
subwavelength apertures, fibers, dipols or atoms). Due to 
Babinct's principle, the model predicts the enhancement 
or supression also in the reflection spectra. According to 
our model, the power (energy) of the light emitted by the 
subwavelength-dimension ensemble of light sources scales 
as the number of light-sources squared, regardless of pe- 
riodicity of the array of sources. Such an effect is not 
unknown one in the physics. The famous Dicke quantum 
model of the superradiance emission of a subwavelength 
ensemble of atoms predicts the same scaling behav- 
ior. Therefore, the mechanism described in the present 
paper could be interpreted as a non-quantum analog of 
the superradiance emission of a subwavelength ensemble 
of coherent light-sources. The evident resemblance be- 
tween our model and the Dicke model also indicates that 
the interference of waves at the detector could lead to 
a new effect, namely the enhancements and resonances 
(versus period of the array) in the total power emitted 
by the periodic array of quantum oscillators (atoms). A 
quantum reformulation of our model, which will be pre- 
sented in the next paper, could also help us to under- 
stand better why a quantum entangled state is preserved 
on passage through a hole array [511 ]. 

In conclusion, we have demonstrated a new mecha- 
nism that could provide a great resonant and nonreso- 
nant transmission enhancement of the light or de Broglie 
particle waves passed through the subwavelength width 
slits not by the surface waves, but by the constructive 
interference of diffracted waves (beams generated by the 
apertures) at the detector placed in the far-field zone. 
The model shows that the beams generated by multiple, 
subwavelength wide slits can have similar phases and can 
add coherently. If the spacing of the slits smaller than the 
optical wavelength, then the phases of the multiple beams 
are nearly the same and beams add coherently (the light 
power and energy scales as the number of light-sources 
squared, regardless of periodicity). If the spacing is 
larger, then the addition is not so efficient, but still leads 
to enhancements and resonances (versus wavelength) in 
the total power transmitted. In contrast to other models, 
the mechanism depends neither on the nature (light or 
matter) of the beams (continuous waves or pulses) nor on 
material and shape of the multiple-beam source (arrays 
of 1-D and 2-D subwavelength apertures, fibers, dipoles 
or atoms). The verification of the results by comparison 
with data published in the literature supports the model 
predictions. The Wood anomalies in transmission spec- 
tra of gratings, a long standing problem in optics, follow 
naturally from the interference properties of our model. 
The new point is the prediction of the Wood anomaly in a 
classical Young-type two-source system. The new mech- 
anism could be interpreted as a non-quantum analog of 
the superradiance emission of a subwavelength ensem- 
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ble of atoms (the light power and energy scales as the 
number of light-sources squared, regardless of periodic- 
ity) predicted by the well-known Dicke quantum model. 
We stress again that the plasmons or trapped electro- 
magnetic modes do not affect the principle of the en- 
hancement based on the classic or quantum interference 
of diffracted waves (beams generated by the independent 
subwavelcngth sources) at the detector placed in the far- 
field zone. The plasmon-polaritons could provide just ad- 
ditional enhancement by increasing the power (energy) 
of each beam. The analytical formulas derived in the 
present study could be useful for experimentalists who 
develop nanodevices based on transmission and beaming 
of light or matter waves by subwavelengths apertures. 

This study was supported in part by the Framework 
for European Cooperation in the field of Scientific and 
Technical Research (COST, Contract No MP0601). 
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